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Abstract 
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1 Introduction 

Molecular electronics and transport have attracted 
much more attention since molecules constitute 
promising building blocks for future generation 
of electronic devices. The transport through 
molecules was first studied theoretically in 1974 [1]. 
Later several numerous experiments [5J O [71 [HI IH] 
have been performed through molecules placed be- 
tween two electrodes with few nanometer sepa- 
ration. Full quantum mechanical treatment is 
required to characterize the transport through 
molecules. Some ab initio calculations [H [31 [3] are 
also used to study the current-voltage characteris- 
tics of a molecule. Electrical conduction through 
isomer molecules strongly depends on the relative 
position of the atoms in molecules and their cou- 
pling strength to the two leads, namely, left lead 
and right lead. Here we describe theoretically the 
electron transport in some specific models of isomer 
molecules. Based on the scanning probe technique 
measurement, conductance of molecular systems is 
directly possible [IHl [IH [H [H [TH [15] . Theoreti- 
cally there exist several formulations P^Bl [T71 [TBI fTUl 
[20l [2T1 [22] for the calculation of conductance based 
on Landauer formula and the seminal 1974 paper 
of Aviram and Ratner [1] . At much low tempera- 
tures and low bias voltages the electron transport 
becomes coherent through the molecule. Here we 
assume that the dissipation and equilibration pro- 
cesses occur only in the two contacting leads and 
this approximation enables to describe the propa- 
gation of an electron by means of single particle 
Green's function. This theory is much more flex- 
ible than any other theoretical approach and also 
applicable to any system described by a Hamilto- 
nian with a localized orbital basis. By using this 
method the electronic transport properties of any 
system can be studied very easily with a very small 
computational cost. In that case we have to know 
only the Hamiltonian matrix for the molecule but 
no need to know anything about the electronic wave 
function. 

Here, we investigate the conductance (g) and 
current- voltage {I-V) characteristics of three differ- 
ent isomer molecules and compute the effect of these 
isomers on the electrical conduction. 

The present article is organized as follows. In Sec- 
tion 2, we describe the formulation of conductance 
g by calculating the transmission probability T and 
current / for any finite size conducting system at- 
tached to two semi-infinite metallic leads by the use 
of Green's function technique. Section 3 describes 



the conductance and current-voltage characteristics 
of single isomer molecules. Finally, we conclude our 
results in Section 4. 

2 Formulation of ^f, T and /: 
Green's function technique 

Here we give a brief description for the calcula- 
tion of transmission probability (T), conductance 
(g) and current (/) through a finite size conducting 
system attached to two semi-infinite metallic leads 
by the use of Green's function technique. 

Let us first consider a one-dimensional conductor 
with N number of atomic sites (array of filled cir- 
cles) connected to two semi-infinite leads, left lead 
and right lead as shown in Fig. [H The conducting 
system in between the two leads can be an array 

Leadl Lead2 

Figure 1: Schematic view of a one-dimensional con- 
ductor with N number of sites (filled circles) at- 
tached to two leads. The first and the last sites are 
labeled by 1 and N, respectively. 

of some quantum dots, or a single molecule, or an 
array of some molecules, etc. At low voltages and 
low temperatures, the conductance of a conductor 
can be written through the Landauer conductance 
formula as [53], 



where g is the conductance and T is the transmis- 
sion probability of an electron through the conduc- 
tor. This transmission probability can be expressed 
in terms of the Green's function of the conductor 
and the coupling of the conductor to the two leads 
in the following way, 

T = Tr [FiG^FflG^] (2) 

where and are respectively the retarded and 
advanced Green's function of the conductor. F^ 
and Tn are the coupling terms between the conduc- 
tor and the two leads. For the complete system, i.e., 
the system with the conductor and the two leads, 
the Green's function is defined as, 

G = {e-Hy' (3) 
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where e ~ E + it]. E is the energy and 77 is a very 
small number which can be put as zero in the lim- 
iting approximation. The above Green's function 
corresponds to the inversion of an infinite matrix 
which consists of the full system. It can be parti- 
tioned into different sub-matrices that correspond 
to individual sub-systems. 

The Green's function for the conductor can be 
effectively expressed in the form, 



Gc — {€ — He — — 'Eji) 



(4) 



where He is the Hamiltonian of the conductor sand- 
wiched between the two leads. The Hamiltonian of 
the conductor in the tight-binding framework can 
be written within the non-interacting picture in this 
form, 



where c| (c^) is the creation (annihilation) opera- 
tor of an electron at site i, e^'s are the site energies 
and t is the nearest-neighbor hopping strength. In 
Eq. HI Sl = h'l^gLhLc and T^r = hRcQRh}^^ are 
the self-energy terms due to the semi-infinite leads. 
gL and gR are respectively the Green's function for 
the left and right leads. Hlc and Hrc are the cou- 
pling matrices and they will be non-zero only for 
the adjacent points in the conductor, 1 and N as 
shown in Fig. [H and the leads respectively. The 
coupling terms T l and T r for the conductor can be 
calculated through the expression [53] , 



{L,R} 



y.r 



(6) 



where F^^ and F|^ j^-j are the retarded and ad- 
vanced self-energies respectively and they are con- 
jugate with each other. Datta et al. [24] have shown 
that the self-energies can be expressed in this form. 
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where A^^^^j are the real parts of the self-energies 
for the two leads. The imaginary parts A^j^ Ry of 
the self-energies represent the broadening of the en- 
ergy level due to the coupling of the conducting 
system with the two leads respectively. By calcu- 
lating some algebra these real and imaginary parts 
of the self-energies can also be determined in terms 
of coupling strength (r^i^ Ry) between the conduc- 
tor and the two leads, injecting energy (E) of the 
transmitting electrons and hopping strength (t) be- 
tween nearest-neighbor sites in the leads. Thus the 



coupling terms F^ and Tr can be written in terms 
of the retarded self-energy as. 



(8) 



Thus by calculating the self-energies due to the left 
and right leads the coupling terms F^ and Tr can 
be easily obtained and then the transmission prob- 
ability (T) will be calculated from the expression as 
mentioned in Eq. [2j 

Since the coupling matrices h^c and Hrc are non- 
zero only for the adjacent points in the conductor, 
1 and N as shown in Fig. (TJ the transmission prob- 
ability can be expressed in the following way, 

T{E, V) - 4Af,{E, V)A'^j^{E, V)\Gin{E, V)\^ 

(9) 

For the sake of simplicity, here we assume that the 
entire voltage (V) is dropped across the conductor- 
electrode interfaces. 

The current through the bridge system can be 
written in the following form [23] . 



EF+eV/2 



T{E, V)dE 



(10) 



eV/2 



where Ep is the Fermi energy of the conductor. Us- 
ing the expression of T{E, V) as in Eq. [5] the final 
form of I{V) will be, 

Ep+eV/2 

m = ^ / Ai,{E,V)A^^{E,V) 



EF-eV/2 

X \GiN{E,V)\^dE 



(11) 



Eq. [U Eq. [5] and Eq. [TT] are the final working 
formulae for the calculation of conductance and 
current-voltage characteristics respectively for any 
finite size conducting system connected to two semi- 
infinite leads. 

Here, we shall describe conductance-energy and 
current-voltage characteristics by using the above 
formulations for some specific models of single iso- 
mer molecules. Throughout the paper, we use the 
units where c = h = e = 1. 

3 Conductance and current- 
voltage characteristics of 
isomer molecules 

Electron transport is strongly affected by the rela- 
tive position of the atoms in single isomer molecules. 



3 



and, to describe this effect we consider three dif- 
ferent models of isomer molecule. Figure [2] corre- 
sponds to the schematic geometry of the three dif- 
ferent isomer molecules. They are respectively de- 
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Figure 2: Schematic representation of three dif- 
ferent isomer molecules: pyridazine, pyrazine and 
pyrimidine those are attached to two leads. 

fined as pyridazine, pyrazine and pyrimidine, where 
relative position of the two nitrogen atoms changes 
accordingly as shown in the figure. These single 
molecules are connected to two semi-infinite leads 
by thiol (S-H) groups (not shown here in the figure) . 
In actual experimental setup, two leads made from 
gold (Au) are used and molecule attached to the 
leads by thiol (S-H) groups in the chemisorption 
technique where hydrogen (H) atoms remove and 
sulfur (S) atoms reside. Here we assume that the 
molecules are connected to the two leads via sulfur 
atoms by removing the extreme left and right hy- 
drogen atoms of each molecule as shown in Fig. [51 
It should be noted that the electron transport 
is strongly affected by the molecule to lead cou- 



pling strength and in this article we shall describe 
the isomeric effect on electron transport both for 
the weak and strong molecule-lead coupling lim- 
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Figure 3: Conductance g as a function of inject- 
ing electron energy E for the weak coupling limit of 
three isomer molecules, where (a), (b) and (c) cor- 
respond to the pyridazine, pyrazine and pyrimidine 
molecules, respectively. 

its. The weak coupling limit is mentioned through 
the condition T^L^ny << t, while the strong cou- 
pling limit is denoted as T[l,r} ^ t- '''{l,r} is the 
molecule to lead coupling strength and t is the hop- 
ping strength between nearest-neighbor site in the 
molecule. Throughout the paper, we set the values 
of these parameters in the weak coupling limit as 
tl = tr = 0.5, t = 3, and, for the strong coupling 
limit these parameters are taken as t^, = = 2, 
i = 3. Now, we try to characterize the conductance- 
energy and current-voltage characteristics of these 
single isomer molecules. 

In Fig. [3l we plot the conductance variation as a 
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function of energy E of these three isomer molecules 
for the weak molecule to lead coupling limit. Fig- 
uresEfa) , (b) and (c) correspond to the conductance 
variation of pyridazine, pyrazine and pyrimidine 
molecules respectively. The conductance shows os- 
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Figure 4: Conductance g as a function of injecting 
electron energy E for the strong coupling limit of 
three isomer molecules, where (a), (b) and (c) cor- 
respond to the pyridazine, pyrazine and pyrimidine 
molecules, respectively. 

dilatory behavior with sharp resonant peaks for 
some particular energy values, while for all other 
energy values it almost vanishes. This behavior 
can be understood as follows. Transmission of an 
electron through a molecule takes place only when 
the incident energy matches with anyone of the en- 
ergy eigenvalues of the molecule. For this particular 
energy value the electron transmits almost ballisti- 
cally from the source to drain, and accordingly, con- 
ductance shows a sharp peak. Now the electrons are 
carried from the left lead to right lead through iso- 
mer molecules and hence the electron waves prop- 



agating along the two arms of the ring may suffer 
a phase shift between themselves, according to the 
result of quantum interference between the various 
pathways that the electron can take. Therefore, 
the probability amplitude of the electron across the 
ring becomes strengthened or weakened. It empha- 
sizes itself especially as transmittance cancellations, 
some peaks do not reach to unity anymore, and anti- 
resonances in the transmission (conductance) spec- 
trum. From Fig. [S^a) and Fig. ^h) it is observed 
that at resonances the conductance shows the value 
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Figure 5: Current / as a function of the applied bias 
voltage V for the weak coupling limit of three iso- 
mer molecules, where the solid, dotted and dashed 
lines are respectively for the pyridazine, pyrazine 
and pyrimidine molecules. 

2, i.e., for these resonances transmission probabil- 
ity T goes to unity (since from Landauer formula 
we get g — 2T hy putting e = h = 1). On the other 
hand. Fig. ^c) shows that the conductance peak 
does not reach the value 2 anymore. Though all 
these three molecules are attached symmetrically to 
the two leads, yet the probability amplitude across 
the pyrimidine molecule decreases compared to the 
other two molecules. So it can be emphasized that 
the electron transmission is strongly affected by the 
relative position of the two nitrogen atoms in these 
isomer molecules. Thus, different molecular struc- 
ture can significantly affect the transport property 
even if they are isomer. 

Another important feature observed from this 
Fig. [3] is that, for the pyrazine (Fig. [3Kb)) and 
pyrimidine (Fig.ISjc)) molecules electron transmis- 
sion starts at very low energy value, i.e., for low 
applied bias voltage V. On the other hand, the 
pyridazine molecule starts electron transmission at 
quite higher energy value (see Fig.[3Ka)). Thus, the 
threshold bias voltage for electron transmission also 
depends on the relative position of the two nitrogen 
atoms in these isomers. 
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To investigate the effect of conductance beliavior 
in strong molecule to lead coupling limit, we plot 
the results of 5 as a function of energy E for these 
three isomer molecules in Fig. |4l where (a), (b) 
and (c) are respectively for the pyridazine, pyrazine 
and pyrimidine molecules. From these curves it is 
observed that the resonant peaks get substantial 
widths and electron transmission takes place for 
wide range of energy values through these molec- 
ular bridge systems. This is due to the substantial 
broadening of the molecular energy levels caused by 
the strong coupling of the molecule to the two semi- 
infinite leads [23] . These effects have strong depen- 
dence on the current-voltage characteristics which 
we are going to discuss in the following parts. 

Now, we study the current-voltage characteris- 
tics of these molecular bridge systems. The cur- 
rent / is computed by the integration procedure 
of the transmission function T. The behavior of 
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Figure 6: Current / as a function of the applied bias 
voltage V for the strong coupling limit of three iso- 
mer molecules, where the solid, dotted and dashed 
lines are respectively for the pyridazine, pyrazine 
and pyrimidine molecules. 

the transmission function T with energy E is ex- 
actly similar (differ only in magnitude by the fac- 
tor 2 due to the existence of the relation g = 2T) 
to that as presented in Fig. [3] and Fig. ID Fig- 
ure [5] represents some current- voltage characteris- 
tics of the three molecules in the weak molecular 
coupling limit, where the solid, dotted and dashed 
lines give the current for the pyridazine, pyrazine 
and pyrimidine molecules, respectively. The cur- 
rent shows staircase-like behavior with sharp steps. 
This is due to the discreteness of the molecular res- 
onances as shown in Fig. [3l From this figure it is 
clearly observed that, both the pyrazine and pyrim- 
idine molecules show finite non-zero value of current 
(dotted and dashed curves) for very small applied 
bias voltages, while the pyridazine molecules gives 



non-zero value of current (solid curve) at higher bias 
voltages. 

In the presence of strong molecule to lead cou- 
pling limit, the staircase-like behavior disappears 
and current varies quite continuously with the ap- 
plied bias voltage. The variation of the current for 
these three isomer molecules in the strong coupling 
hmit is plotted in Fig. [SI where the soHd, dotted and 
dashed lines correspond to the same meaning as ear- 
lier. For this coupling limit, molecular resonances 
get substantial width (Fig. [4]) and since the cur- 
rent is computed by the integration procedure of the 
transmission function T, it gets a continuous varia- 
tion with the applied bias voltage V. Similar to the 
weak coupling case, the pyrimidine molecule gives 
non-zero value of the current at sufficient higher 
bias voltages compared to the other two molecules 
in this strong coupling case. Another important ob- 
servation is that, in the strong coupling limit cur- 
rent amplitude becomes very large than that of the 
current amplitude in the weak coupling limit for 
these molecular bridge systems. 

4 Concluding remarks 

To conclude, we have investigated the effect of iso- 
mers and molecule to lead coupling strength on 
electron transport through some models of different 
types of isomer molecules. Depending on the rela- 
tive position of the two nitrogen atoms in these iso- 
mer molecules, threshold bias voltage for electronic 
conduction changes which provides an important 
signature for the fabrication of different molecular 
bridge systems. The conductance shows fine reso- 
nant peaks which gives the staircase-like behavior 
with sharp steps in current for the weak coupling 
limit, while in the strong coupling limit, the reso- 
nant peaks get substantial widths, and accordingly, 
the current varies almost continuously with the ap- 
plied bias voltage. 
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